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ABSTRACT 


A wind tunnel investigation of low speed flow over flapping airfoils and airfoil combi¬ 
nations was performed using flow visualization and laser Doppler velocinretry. Specifically, 
three cases were studied: A NACA0014 airfoil oscillating in a sinusoidal plunge mode, A 
NACA0014 airfoil oscillating in a sinusoidal plunge mode near a ground plane, and two 
NACA0014 airfoils arranged in a biplane configuration and oscillating in counterphase in a 
sinusoidal plunge mode. The plunge amplitude-to-airfoil chord ratio was 0.4, the reduced 
frequency of oscillation was 1.0 and the Reynolds number based on airfoil chord was set at 
8760. 

Conditionally sampled measurements of the axial flow velocity were taken at nu¬ 
merous flow field points providing detailed information about the flow features generated 
by this type of flapping motion. These measurements were complemented by tinre-averaged 
flow field data and by visualization of the instantaneous flow field at various points during 
the flapping cycle. Furthermore, the thrust generated by the sinusoidal plunge motion was 
measured with a laser range finder. The results shows that vortex shedding occurs both 
from the airfoil leading and trailing edge. 
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I. INTRODUCTION 


A. BRIEF HISTORY OF FLAPPING WING RESEARCH 

The desire “to fly like the birds” is almost as old as humanity itself. Indeed, Leonardo 
da Vinci was so fascinated by the flapping flight phenomenon that he made many sketches 
of bats and summarized his flapping wing studies in a book manuscript. However, little 
progress was made during the following centuries until the German flight pioneer Otto 
Lilienthal undertook systematic glider flight tests with flapping wings and published his 
findings in 1889. Knoller [17] and Betz [2] appear to have been the first ones to offer 
an explanation for the development of thrust due to wing flapping. Katzmayr [15] first 
presented wind tunnel experiments which showed that a wing mounted in an oscillating air 
flow experienced a thrust force. Birnbaum [3] first succeeded to develop a theory which 
correctly predicted the lift and thrust forces developed by flapping wings. In the following 
decades Kuessner and Theodorsen [27] were able to eliminate the restriction of Birnbaum’s 
theory from low-frequency oscillations. This Theodorsen-Kuessner oscillatory thin-airfoil 
theory thereafter became the standard tool to analyze airfoil flutter problems. Garrick [8], 
applied Theodorsen’s theory to the analysis of flapping wing propulsion, which remains the 
classical reference work to this day. 

The rapid development of fixed-wing airplanes in conjunction with the obvious me¬ 
chanical complications introduced by flapping wings soon discouraged the further develop¬ 
ment of ornithopters. Indeed, the demonstration of a successful manned ornithopter is still 
an aeronautical engineering challenge which is being pursued by flapping wing enthusiasts, 


1 



such as DeLaurier [6] and his team at the University of Toronto, among others. Further 
discouragement came from Garrick’s results which showed that flapping wings have low 
propulsion efficiencies unless oscillated at very low frequencies. However, Schmidt [25] tried 
to overcome this concern during and after World War II by developing a device which he 
called ’’wave propeller”. He argued that the vortices generated by a flapping wing carry 
energy which can be converted into additional thrust by a wing mounted in the wake of a 
flapping wing, The aft-wing thus is exposed to an oscillatory flow which generates thrust by 
virtue of the Katzmayr effect. Schmidt claimed to have achieved wave propeller efficiencies 
comparable to conventional propeller efficiencies. 

Nevertheless, the practical development potential of flapping wings for either pure 
propulsion or as an integrated lift/propulsion system was regarded as unattractive until very 
recently. Flapping wing studies therefore largely remained restricted to scientists interested 
in bird flight or fish propulsion problems. 

An unexpected revival of interest in the systematic study of flapping wing flight 
phenomena occurred only a few years ago with the announcement of a major initiative by 
the Defense Advanced Research Projects Agency (DARPA) to encourage the development 
of micro-air vehicles (MAVs). The use of flapping wings for vehicles with dimensions not 
exceeding six inches in length or span is an obvious option because of the low efficiency of 
conventional propellers at low Reynolds number. Therefore, this DARPA initiative sparked 
a number of investigations which have been reviewed in June 2000 at a special symposium 
held at the University of Notre Dame to be published by AIAA in 2001. This publication 
is referred to for a detailed review of the current status of these studies and we therefore 
limit ourselves to a brief description of the major developments in recent years. 

Of special note is the replacement of the Theodorsen-Kuessner oscillatory thin-airfoil 
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theory by an approach which enables the computation of incompressible potential flow past 


oscillating airfoils of arbitrary shape. This is being accomplished by the placement of sources 
and vortices on the airfoil surface rather than along the chord line. This so-called panel 
method was pioneered by Giesing [9], who generalized the method of Hess and Smith [10], for 
steady airfoil flow. It was further developed at the Naval Postgraduate School, as described 
by Teng [26], Platzer et al [22], and Jones et al [16] . The rapid increase in computer power 
over the past few years also made it possible to apply Navier-Stokes computations to the 
analysis of flapping airfoil aerodynamics, thus making it possible to compute strong viscous 
effects, as shown by Tuncer and Platzer [30], Tuncer et al [29], Tuncer and Platzer [28], Isogai 
et al [11], and Ramamurti and Sandberg [24]. Also, a number of investigators succeeded to 
analyze three-dimensional effects. Of special note here is the early work of Lan [19], who 
applied the vortex lattice method and showed that the interaction between a flapping fore 
and aft-wing can be quite beneficial. Peng-fei Liu [20] developed a three-dimensional panel 
code and, most recently, Neef (private communications) even succeeded to analyze three- 
dimensional flapping wing effects using a Navier-Stokes code. Experimentally, Freymuth 
provided valuable flow visualizations of the vortical wakes generated by plunging or pitching 
airfoils. Jones et al [12] and Lai and Platzer [18] showed in their flow visualization studies 
that the vortical wake becomes asymmetric as soon as the reduced plunge velocity exceeds 
a critical value. Jones et al [16], and Anderson et al [1] made thrust measurements which 
provide information about the dependence on flapping frequency and amplitude. 
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B. THE PHYSICS OF THRUST GENERATION DUE TO 


WING FLAPPING 

As already mentioned, Knoller [17] and Betz [2] were the first ones to offer an 
explanation for the birds’ ability to generate a propulsive force by means of flapping their 
wings. Consider the airfoil in sinusoidal flapping while also flying forward. As the airfoil 
moves through its mean position during the downward stroke, it is effectively exposed to 
a flow with positive angle of incidence. Similarly, it sees a negative incidence angle during 
the upstroke. If, for simplicity, the resulting aerodynamic force is assumed to be essentially 
perpendicular to the instantaneous approach flow angle, then decomposition into a force 
component parallel to the flight velocity vector will produce a small sinusoidally varying 
thrust force. 

It is understood that this explanation is greatly simplified. The actual flow which is 
produced is considerably more complicated, even if only attached inviscid flow is assumed. 
As is well known, every change in incidence will produce a starting vortex which is shed from 
the trailing edge. Therefore, a sinusoidally oscillating airfoil will generate a vortex street 
behind the airfoil. This phenomenon can be simulated quite easily with the above mentioned 
unsteady panel codes. The vortex street consists of an upper row of counterclockwise 
vortices and a row of lower clockwise vortices Fig. 1.1. This vortex street therefore is 
just the opposite of the well known Karrnan vortex street. If time-averaging is applied at 
some location cutting the wake in the normal direction to the free-stream a jet profile is 
obtained, whereas the same operation applied to the Karrnan vortex street would produce 
a wake profile, Fig. 1.2 . The vortex street produced by the flapping foil in effect produces 
a jet flow. This is to be expected since the thrust experienced by the airfoil must be found 


4 




Figure 1.1. Thrust-Indicative Wake Pattern, after [12] 


as momentum increase in the fluid. As shown by Jones et al [12], this jet flow can indeed be 
measured and is in good agreement with the panel code predictions. Jones et al also found, 
however, that the symmetric vortex street switches into a deflected street as soon as the 
non-dimensional plunge velocity exceeds a critical value, again in good agreement between 
the flow visualizations and the panel code predictions. 

This encouraging agreement between the measurements and the inviscid flow pre¬ 
dictions of Jones et al [16] might give the impression that the physics of flapping airfoils is 
understood reasonably well and that the prediction of the achievable thrust can be made 
with considerable confidence. Unfortunately, even disregarding the three-dimensional flow 
effects introduced by finite-span wings, the range of validity of inviscid flow predictions is 
severely limited by the onset of dynamic stall. This seems to be particularly true at the low 
Reynolds numbers typically required for micro-air vehicles. 
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Figure 1.2. Drag-Indicative Wake Pattern, after [12] 


C. THE PHENOMENON OF DYNAMIC STALL 

Helicopter aerodynamicists are quite familiar with the phenomenon of dynamic stall. 
A helicopter blade in forward flight can experience dynamic stall while it is in the so-called 
retreating blade position, exposing the blade to high incidence angles. For this reason, 
the aerodynamics of pitching airfoils experiencing dynamic stall has been studied to a 
considerable extent both experimentally and computationally. The current state-of-the-art 
can be found in the recent review papers of Carr and Chandrasekhara [4], and Ekaterinaris 
and Platzer [7]. During the downstroke, a stall vortex forms near the leading edge which 
moves over the upper surface and eventually leaves the airfoil, strongly interacting with the 
trailing edge vortex which forms at the trailing edge. As soon as the stall vortex leaves the 
airfoil, rapid lift loss and changes in pitching moment occur, leading to a severe hysteresis 
loop in lift, drag and pitching moment. 

Most dynamic stall studies have been limited to pitching airfoils because of the 
importance of this motion for helicopter blades. However, for thrust generation it is well 
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known (and has been shown in detail by Jones and Platzer [14]) that a pure plunging or 
a combined pitch/plunge motion is required in order to produce significant thrust forces. 
Very few experiments involving dynamic stall due to pure plunge have as yet been carried 
out, but several Navier-Stokes computations have been completed by Tuncer and Platzer, 
Isogai et al, and Ramamurti and Sandberg, as already noted above. These computations 
clearly show the possibility of dynamic stall for values of flapping frequency, amplitude, and 
Reynolds number typical for MAV flight. 

D. MOTIVATION FOR THIS INVESTIGATION 

Because of the very limited experimental information available for plunging airfoils 
it is the objective of this investigation to conduct detailed experiments which provide flow 
visualization and quantitative measurements of the flow near and downstream of flapping 
airfoils. As shown by Jones and Platzer [14], the performance of flapping airfoils is enhanced 
if two airfoils are flapping in counterphase, effectively producing the favorable ground effect 
experienced by birds flying near the water surface. Therefore, three experiments are carried 
out, 

1. Two airfoils plunging sinusoidally in counterphase. 

2. A single airfoil plunging sinusoidally near a ground plane. 

3. A single airfoil plunging sinusoidally. 

so that the differences caused by these three configurations can be understood in 
sufficient detail. 
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II. COMPUTATIONAL ANALYSIS 


A. NAVIER-STOKES ANALYSIS 

In the companion computational study by Castro [5], a thin-layer Navier-Stokes 
code was used to compute the two dimensional flow over flapping airfoils. Two cases of his 
study were investigated in the present work, with Laser Doppler Velocimetery (LDV) and 
flow visualization. These two cases are: 

• Biplane configuration in pure plunge motion, k=l, = 0.1, Re = 10000, a = 0°, 
h=0.4 

• Single wing configuration in pure plunge motion, k=l, M ,^ = 0.1, Re = 10000, a = 0°, 
h=0.4 

His results are introduced in the next two sections, mainly to provide some insight 
into the basic flow features which one expects to find in the experiments. However, it is 
important to recall that the computations are based on a fully laminar flow assumption. 

B. BIPLANE CONFIGURATION 

The wings’ flapping motion is presented in Fig. 2.1. In the biplane case, the flow 
behavior around the upper airfoil is the same as around the lower one. It is antisymmetric 
about the plane of symmetry. Hence, from now on, the flow around the upper airfoil will 
be the one that will be described and the word airfoil will mean the upper one. 

The non-dimensional velocity around the biplane configuration is presented as vec- 
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tors at 8 different plunge phase angles 4>, namely 0°, 45°, 90°, 135°, 180°, 225°, 270°, and 


315°. The non-dimensional velocity is presented as vectors, while the vorticity and the 
pressure distribution are presented as flood contours in Fig. 2.2, Fig. 2.3, Fig. 2.4, Fig. 2.5, 
Fig. 2.6, Fig. 2.7, Fig. 2.8, and Fig. 2.9. 

The flapping cycle starts at = 0°, Fig. 2.2. The wing is in the uppermost position 
and starts to move toward the plane of symmetry. At this position, = 0°, the flow near 
the lower surface of the airfoil has two counter-clockwise vortices surrounding a clockwise 
vortex. These vortices were generated during the upstroke, as will be explained in more 
detail. The velocity near the upper surface is close to the freestream velocity. In the near 
wake region (half period, since k=l, c=l, hence L = 2i r), the velocity decreases. Further 
downstream there is a clockwise vortex which originated from the dynamic stall vortex shed 
during the previous cycle. 

As the airfoil moves down to = 45°, the effective angle of attack starts to increase 
on the upper surface causing a small pressure drop, Fig. 2.3. Also, the stagnation point 
moves to the lower surface close to the leading edge. Meanwhile, the three vortices which 
were sitting near the lower surface start to get shed downstream, such that the first one 
moves to the near wake region causing the velocity to increase just behind the trailing edge. 
The pressure increases between the two airfoils as the velocity decreases. Meanwhile, in the 
far wake region the vortex keeps moving downstream. 

The airfoil reaches the mean flapping position at = 90°, Fig. 2.4. By then all the 
vortices are completely shed downstream away from the airfoil surface. The effective angle 
of attack increases and the stagnation point moves further downstream at the lower surface. 
The velocity increases on the upper surface within the first quarter chord and a vortex starts 
to be formed. The pressure increases between the two airfoils causing a velocity increase 
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just behind the trailing edge which helps to shed the vortices into the near wake region. 


In the far wake, the velocity starts to decrease below the free stream value and a region of 
higher pressure starts to be generated near the plane of symmetry. 

As the airfoil keeps moving down, <f> = 135°, the vortex generated on the upper 
surface becomes stronger, Fig. 2.5. As a matter of fact, four vortices start to develop, two 
clockwise and one counter-clockwise within the first half chord. A smaller clockwise vortex 
is generated close to the trailing edge. The velocity decreases on the lower surface and the 
pressure increases to cover an area which extends to half chord behind the trailing edge. 

As the airfoil reaches the lowest position in the flapping stroke, $ = 180°, Fig. 2.6, 
the four vortices on the upper surface become fully developed and the airfoil is in complete 
dynamic stall. The velocity between the two airfoils is slightly higher than the free stream 
value. In both the near and far wake regions, two vortices can still clearly be seen. The 
stagnation point has moved forward toward the leading edge but is still on the lower surface. 
The pressure is a little lower than the free stream value. 

The down stroke ends at $ = 180°. 

As the airfoil starts to move up, the effective angle of attack switches from positive 
to negative values. Hence, the stagnation point starts to move to the upper surface which 
can be seen at = 225°, Fig. 2.7. The vortices on the upper surface start to be shed 
downstream clearing the trailing edge. The velocity increases in the area between the two 
airfoils causing a drop in pressure. In the near wake region, the shed clockwise vortex 
behind the trailing edge causes a drop in velocity near the plane of symmetry. The velocity 
starts to increase gradually toward the far wake. 

As the airfoil approaches the mean position during the upstroke, $ = 270°, Fig. 2.8, 
the effective angle of attack increases causing the velocity to increase near the lower surface. 
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A shed clockwise vortex is seen in the near wake which creates a low pressure zone that 


draws the flow toward it. Meanwhile, this vortex slows the flow velocity dramatically near 
the plane of symmetry. As the flow approaches the far wake, the velocity increases slightly 
near the plane of symmetry, but away from the plane of symmetry the flow direction changes 
in a wavy form due to the existence of the vortex at one chord behind the airfoil and the 
shed one at the 4.5 chord length. 

As the airfoil moves upward during the upstroke, the effective angle of attack keeps 
increasing. This becomes obvious at = 315°, Fig. 2.9. As a result, the velocity increases 
on the lower surface which creates three counter-clockwise vortices and a small clockwise 
one. The stagnation point starts to move downstream on the upper surface. The pressure 
decreases a little in the area between the two airfoils. The lower surface of the airfoil is 
almost in complete dynamic stall. 

As the airfoil arrives at the uppermost point,$ = 360° = 0°, Fig. 2.2, the flow near 
the lower surface of the airfoil is in complete dynamic stall. Multiple vortices exist on the 
lower surface while the fourth is already shed downstream. 

C. SINGLE FLAPPING WING IN PURE PLUNGE 

Similarly, the flow around a single flapping wing in pure plunge was computed. 
The flow fields are shown in Fig. 2.10, Fig. 2.11, Fig. 2.12, Fig. 2.13, Fig. 2.14, Fig. 2.15, 
Fig. 2.16, and Fig. 2.17 at 8 plunge phase angles, namely, 0°, 45°, 90°, 135°, 180°, 225°, 
270°, and 315°. 

The airfoil flaps in pure sinusoidal plunging motion. It flaps in the same fashion as 
the top airfoil in Fig. 2.1. 


12 




Figure 2.1. The Flapping Motion in the Biplane Configuration 
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Figure 2.2. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 0°. 
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Figure 2.3. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 45°. 
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Figure 2.4. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 90°. 
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Figure 2.5. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 135°. 
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Figure 2.6. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 180°. 
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Figure 2.7. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 225°. 
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Figure 2.8. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0° , Moo = 0.1, k = 1, Reoo = 1 x 10 4 , 4> = 270°. 
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Figure 2.9. Vorticity and Pressure Distribution for Biplane Configuration, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, Re^ = 1 x 10 4 , $ = 315°. 
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The flapping cycle starts with the wing in the uppermost position at $ = 0° 


Fig. 2.10. At that position, the stagnation point is located near the leading edge. Two 
clockwise and one counterclockwise vortices exist near the lower surface of the airfoil. These 
vortices were generated during the previous upstroke. 

As the airfoil starts the downstroke, $ = 45°, the effective angle of attack increases 
causing the stagnation point to move to the lower surface behind the leading edge, Fig. 2.11. 
The velocity increases on the upper surface causing a slight decrease in pressure. The three 
vortices which existed on the lower surface have moved downstream and two of them have 
already cleared the airfoil surface. The pressure on the upper surface is less than that on 
the lower surface. 

As the airfoil approaches the mid-flapping position at $ = 90°, the effective angle of 
attack keeps increasing and the stagnation point moves further downstream from the leading 
edge on the lower surface, Fig. 2.12. On the upper surface a vortex starts to develop within 
the first quarter chord. Meanwhile all vortices have been shed away from the lower surface. 
The pressure value keeps increasing on the lower surface especially within the quarter chord 
behind the leading edge. 

As the airfoil keeps moving down during the downstroke, $ = 135°, the vortices on 
the upper surface keep growing, Fig. 2.13. Two clockwise vortices and one counter-clockwise 
vortex lying between them can be well distinguished. Another vortex is generated on the 
upper surface just before the trailing edge. 

As the airfoil approaches the lowest point in the downstroke, $ = 180°, the flow near 
the upper surface is in full dynamic stall, Fig. 2.14. The three vortices on the upper surface 
are well developed while the one that was at the trailing edge has been shed downstream. 
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The stagnation point is located near the leading edge. The pressure on the upper surface 


is less than on the lower surface. 

As the wing starts the upstroke, the effective angle of attack changes from positive 
to negative and the stagnation points starts to move downstream on the upper surface, 
$ = 225°, Fig. 2.15. The velocity increases on the lower surface causing a slight drop in 
pressure. Meanwhile, on the upper surface, the pressure starts to build up. The previously 
generated vortices start moving toward and past the trailing edge. 

As the airfoil approaches the mid-flapping position during the upstroke, $ = 270°, 
Fig. 2.16, the negative effective angle of attack keeps increasing causing the stagnation point 
to move downstream on the upper surface. A high pressure zone is well developed on the 
upper surface and concentrated mainly around the stagnation point within the first quarter 
chord. The velocity increases on the lower surface causing a pressure drop and a vortex 
starts to develop. 

As the airfoil keeps moving during the upstroke, $ = 315°, Fig. 2.17, the velocity 
keeps increasing on the lower surface. Four vortices are now completely developed. 

As the wing reaches the uppermost point, it starts a new cycle as presented at 
$ = 0°, Fig. 2.10 
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Figure 2.10. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0° , Moo = 0.1, k = 1, Re = 1 x 10 4 , 4> = 0°. 
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Figure 2.11. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 45°. 
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Figure 2.12. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 90°. 
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Figure 2.13. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 135°. 
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Figure 2.14. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 180°. 
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Figure 2.15. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 225°. 


29 









































































































Figure 2.16. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 270°. 
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Figure 2.17. Vorticity and Pressure Around a Single Wing in Pure Plunge, NACA 0014, 
a = 0°, Moo = 0.1, k = 1, R eoo = 1 x 10 4 , 4> = 315°. 
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III. EXPERIMENTAL TOOLS 


A. GENERAL DESCRIPTION OF THE SYSTEM 

1. Wind Tunnel Description 

The low-speed, continuous, in-draft wind tunnel used in the present investigation is 
shown in Fig. 3.1. Air is ingested from inside the building through a square 4.5 x 4.5 nr 
inlet, converging through a 9:1 bell shaped contraction to a 1.5 x 1.5 nr test section. Tunnel 
speed is controlled by a variable pitch fan driven by a constant speed electric motor. Motor 
and fan vibration are isolated from the test section by rubber sleeves on each side of the 
nrotor/fan assembly. The tunnel velocity range is 0 to 9.5 nr/s. 

2. Pitot Static Tube 

A pitot static tube was installed above the model for general adjustments of tunnel 
velocity. The pitot-static tube is connected to a MKS Baratron type 223B differential 
pressure transducer. The transducer provides a voltage that is linear with the differential 
pressure, yielding 1 volt at 1 torr. The flow-speed is then given by Eq. (3.1) 


U OQ = 



(3.1) 


The transducer has a reported accuracy of 0.5% of the full-scale reading which, due 
to the nature of Eq. (3.1), results in rather large errors in the measurement of low velocities 
but quite small errors in the measurement at high velocities. The pitot-static tube has been 
calibrated using LDV, and was shown to be sufficiently accurate. 
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Figure 3.1. Wind Tunnel Schematic Diagram 


The pitot static tube has errors associated with the measurement of both the static 
and total pressure. The pitot tube used here has an outer diameter of 8mm, 8 static ports 
aligned symmetrically, 64mm (8 diameters) downstream of the tip base, and a stem approx¬ 
imately 120nrm (15 diameters) downstream of the static ports. According to Pope [23] the 
geometry of the probe should yield about a 0.5% overprediction of the static pressure. 

3. Laser Doppler Velocimetry (LDV) 

LDV is an optical method that utilizes the Doppler principle, measuring the fre¬ 
quency shift of laser light scattered by a particle in the flow. A dual-beam, on-axis, 
backward-scatter system is used to measure the air velocity and turbulence intensity in 
the NPS low-speed, open circuit smoke tunnel. 

In the dual-beam arrangement utilized, a beam of laser light is split in two, the 
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frequency of one of the beams is shifted, and both beams are transmitted such that they 


cross each other, thereby generating the so-called probe volume. Particles moving through 
the probe volume scatter the laser light. The scattered light is collected, and sent via an 
optical receiving fiber to a Photo Multiplier Tube (PMT). A photo detector in the PMT 
converts the light pulse into a voltage signal at the Doppler frequency. The signal produced 
by the photo detector is sent to the IFA 755 signal processor via a frequency downshifter. 
The signal processor determines the Doppler frequency by performing a burst spectrum 
analysis of the signal. The Doppler frequency is then sent to computer via direct memory 
access (DMA) interface. The PACE software records the signal for later processing into 
velocity information. The LDV probe is mounted on a 3-axis traverse mechanism outside 
the tunnel at the test section window. The LDV probe is fitted with a 750 mm lens allowing 
LDV measurements to the center of the tunnel. 

4. LDV Apparatus 

The LDV system uses a 5W, water-cooled Coherent Innova Series Argon-ion laser 
for the light source. A TSI Model 9201 Colorburst multicolor beam separator unit is used 
to separate the laser colors (green - 514.5 nm, blue- 488 nm, and purple - 476 nm ), split 
the beams and frequency shift one beam in each color. After splitting the beams, one of the 
two blue beams and one of the two green beams were frequency shifted by 40 MHz. The 
purple channel was blocked. The beams were coupled into a 2-conrponent, 4-beam, Model 
9832 fiber optic probe through a Model 9271 coupler (one for each beam). The beams were 
transmitted to the probe via optical receiving fiber. Probe beam spacing is 50 mm and the 
lens focal length is 750 mm. A multi-mode receiving fiber located in the probe received and 
forwarded the reflected light to a Model 9160 photo multiplier tube (PMT). Although the 
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probe is a two-component unit, only one component of the receiving optics was available. 


The output of the PMT is connected to a 9186A frequency downshifter. The downshifter is 
connected to a TSI Model IFA-755 signal processor. The signal processor distinguishes the 
signal burst from the noise based upon the signal-to-noise-ratio. The processor determines 
frequency, phase, burst transit time, and burst arrival time and sends these signals to the 
computer. The DMA board and PACE windows software package on the computer converts 
the signals into velocity, and is used for data display and storage. 

5. Flow Seeding 

Flow seeding was performed using a TSI 9306 six-jet atomizer, operating with com¬ 
pressed air producing particles with approximately 1 micron in size. 

6. Model Description 

For the steady flow measurements, a wing with a NACA 0014 airfoil was fixed to a 
mechanism, which allows variable angles of attack. The wing chord is 64 mm with a span 
of 635 mm. 

For the unsteady flow and thrust measurements, a flapping wing mechanism was 
used. It was designed by Dr. Kevin Jones at the Naval Postgraduate school in 1997. 
Its construction and dynamics are described in [13] from which this description is taken. 
Isometric, profile and top views are shown in Fig. The structure of the fuselage is made of 
aluminum, the front and rear nacelles of balsa-wood, and the wings of wood. The device is 
configured to flap one or two wings. When flapping two wings, the wings are positioned one 
above the other. The two flapping wings are configured to execute a two degree of freedom 
motion, with both pitch and plunge amplitudes adjustable. The phase angle between pitch 
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notch in rear nacelle 



Figure 3.2. Isometric View of Upturned Flapping Mechanism. 


and plunge is also adjustable. Stationary wings may be attached before or after the flapping 
wings in the slots depicted, Fig. 3.2. 

The flapping wings are attached to moving beams, which are driven by eight actuator 
wheels. A pin attached to each of the eight wheels rides in a slot at the end of each beam. 
This attachment changes the circular motion of the wheels into the flapping motion of the 
beams. The beams may be attached to the wheels at a selected wheel radius to control 
plunge amplitude, Fig. 3.3 and Fig. 3.4. Wing pitch is controlled by adjusting the phase 
between the fore and aft wheels. The eight wheels are attached to four worm gears. The 
gears mesh with two worms that are driven by an Astro-Flight Cobalt 40 model airplane 
motor. The motor is powered by a variable volt age/cur rent power supply. Lund [21] 
recommended that the motor voltage should be limited to 26V for a maximum of one 
minute. This voltage limit restricted maximum flapping frequency to 7Hz during the thrust 
measurements. However, for the unsteady flow measurements around the flapping wing, 
longer run times were needed to be able to get enough samples. Therefore, the flapping 
frequency was limited to 5Hz. 

This 26V corresponds to a motor rpm of 21,000 with a 50:1 gear reduction. Current 
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Figure 3.3. Top View of Flapping Mechanism. 



actuator bearinged flapping flapping motor 
wheels pins wings beams 

Figure 3.4. Profile View of Flapping Mechanism. 
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Table 3.1. Flapping Model Parameters 


Wing Configuration 


Single, Bi-wing flapping, no stationary wings 

Wing Type 


Symmetric, approximating NACA 0014 

Wing Span 

b* 

1270 mm 

Fuselage Width 


70 mm 

Useful Lifting Span 

b 

1200 mm 

Wing Chord 

c 

64 mm 

Aspect Ratio 

AR = b 2 /S = b/c 

19.84 full span, 18.75 useful span 

Taper Ratio 

O 

O 

II 

•< 

1 

Sweep Angle 

A 

0 

Mean Wing Separation 


1.4c 

Plunge Amplitude 

h 

h = 0.4 = 25.4 mm 

Pitching Motion 

A a 

None, a geom = 0° 


drawn was always less than 5A. A wing section with a NACA0014 airfoil was used in all the 
experimental work. The wings were manufactured by the Miniature Aircraft Supply model 
helicopter company and were made of balsa wood and laminated hardware, and covered 
with a thin, smooth plastic sheet. The wings can be bolted to the flapping beams, or 
attached as stationary wings in the slots upstream or downstream of the flapping wings. 
The stationary wing feature is useful for experiments to explore Schmidt’s wave propeller. 
The parameters of the flapping mechanism is summarized in Table 3.1 
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B. THE THRUST MEASURING TECHNIQUE 


The experimental arrangement for measuring thrust is depicted in Fig. 3.5. The 
model is suspended with four thin cables 112.8 cm in length. The cables allow displacement 
along the streamwise axis, but hinder motion in other axes. In order to determine thrust, 
the model displacement due to flapping is required. Drag induced by the tunnel flow 
deflects the model rearward. Flapping the wings of the model moves it forward. The net 
displacement, measured by a laser rangefinder, positioned on a two-axis traverse behind the 
model, reflecting the beam from a notch in the model’s rear nacelle, is a measure of the 
thrust. 


After setting tunnel velocity, the laser rangefinder was positioned at its zero voltage 
center point, 13 cm behind the model. Since the model’s rearward displacement due to drag 
increased with velocity, the rangefinder was reset to its center point at each tunnel velocity. 

In order to obtain the desired flapping frequency, the strobe light was set to the 
required frequency and the motor voltage adjusted until the beams flapping the wings 
appeared to be stationary. The model was equipped with a micro switch which triggers 
a square signal once per revolution. This signal was used to obtain the exact flapping 
frequency of the model through the digital oscilloscope. 

During flapping, the laser rangefinder voltage was recorded concurrently with the 
frequency micro-switch voltage by the digital oscilloscope. The rangefinder and frequency 
micro-switch data were recorded with a sample rate of 1000 Hz for a period of approximately 
32 seconds. This ensured collecting enough data. 
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Figure 3.5. Thrust measurement arrangement 


C. UNSTEADY FLOW MEASUREMENTS 


Measuring the instantaneous axial velocity component, u, was essential to know the 


flow characteristics during the wing flapping cycle. The Rotary Machinery Resolver, RMR, 


along with the micro switch circuit allowed conditional sampling for the velocity picked up 


by the LDV system. 


D. TRIGGERING SIGNAL 

A triggering signal was generated by a micro switch which sends a one per revolution 
pulse. An electric circuit was built to get a clean square signal. The RMR uses the rising 


part of the square signal as a reference for position ” ZERO” and the following rising signal 
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as position ’’ZERO”. In the present work, the time interval between the two signals was 


divided into 360 bins which were used to store the velocity measured by the LDV. 

Ideally, the position ’’ZERO” should represent the wing at zero amplitude. But 
due to space limitations inside the model, it was impossible to meet this condition. So, 
the relation between the triggering signal and the actual wing position was required to 
obtain meaningful results. The laser range finder was installed underneath the wing and 
the model set to flap for a while and a set of flapping cycles along with the triggering signal 
was recorded. The obtained results are presented in Fig. 3.6. 

It was found that the trigger signal was generated such that it will trigger at 
(j) = 95.17° considering that the position ZERO is the position when the wing is in the 
uppermost position, which was the convention used in this work. In the present work, the 
maximum amplitude was 0.4C which equals 0.0256 [nr]. Hence, the triggering signal was 
initiated when the wing was at a plunge distance of 0.0023 nr after passing the mean plunge 
position during the downstroke. 

E. CONVENTIONS 

In the present work the following conventions were used: 

• The wing starts the flapping cycle, $ = 0, from the uppermost position. 

• The airfoil leading edge was always placed at X/C = 0, and the trailing edge was 
placed at X/C = 1 

• In the steady measurements, the airfoil chord was placed at Z/C = 0 

• In the unsteady measurements, the mid flapping position was placed at Z/C = 0. 


42 






Figure 3.6. Trigger signal and the plunge motion 
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IV. FLOW VISUALIZATION 


It was of interest to visualize the flow around the NACA0014 airfoil which was used 
in the present work. The flow visualization was done by generating smoke upstream of the 
wing and taking pictures by still and video cameras. A simple mechanism was built to 
change the angle of attack. 

A. SMOKE WIRE SYSTEM 

Flow visualization was performed using a smoke wire, with 0.3 mm diameter Ni- 
Chronre (NI80_CR20) wire stretched between two posts, and spring loaded to account for 
the wire elongation when heated, Fig. 4.1. The resistive wire was heated by applying a 
voltage across its length. The voltage was provided by a Variac set to provide 15-20V. The 
voltage setting was dependent upon tunnel speed and fog fluid flow rate. 

Smoke was generated by dripping Roscoe fog liquid down the wire prior to heating 
it. Supplying the fog fluid was done through a brass tube welded to the fixture at one end of 
the wire. A clear silicon tube was attached to this brass tube which delivered the fog fluid 
from a pressurized container outside the test section. The container was pressurized using 
a bike pump through one of two adapters bolted through the jug’s cap. The transparent 
silicone tube was attached to the second adapter and fog fluid flow was controlled by a pair 
of henrostat clamps. 

To generate streaklines, thin copper wire (0.25 mm diameter) was wrapped around 
the NLChrome wire at intervals and soldered to keep it in place. For the smokewire used, 
the wire/solder beads were spaced at 1 cm intervals, Fig. 4.1. The fog fluid tended to pool 
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on the beads, and it took some time for the wire to burn off these deposits, providing 5 to 


10 seconds of continuous smoke. Close inspection shows that most of the streaklines are 
really two, closely spaced lines. Typically, the fluid pools above and below the beads, and 
it is these two deposits that burn and create the pairs of lines. The nickel chromium wire 
was fixed to a wing with airfoil NACA0014 by two posts such that a 2” distance was kept 
between the wire and the wing to avoid any interference. The flow was illuminated using a 
Smith-Vector Corp 750W Photographic Light. 

A Sony DCR-VX1000 digital video camera was used to capture video images of the 
model in flapping motion and the wing at different angles of attack. A video capturing card 
was used to go through the images and grab the best pictures to present the flow behavior. 
Also, a Kodak DC 260 digital camera was used to capture some of the images. 


B. CASE STUDIES 

Two main study cases were investigated by flow visualization. The first was the 
fixed single NACA0014 airfoil at 0° and 4° angle of attack. The second case was the biplane 
configuration, which was investigated for three different cases, as will be explained later. 

1. NACA0014 Airfoil At Zero Angle Of Attack 

The angle of attack was set to zero, and the tunnel speed to Uoo = 2.2Am/s which 
means that the Reynolds number was in the order of 10,000. The flow around the airfoil 
is presented in Fig. 4.2. A small separation area was noticed near the trailing edge at the 
upper surface. 
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Figure 4.1. Smoke Wire Arrangement With the Copper Nodes on the Smoke Wire 
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2. NACA0014 Airfoil at 4° Angle Of Attack 

At the same tunnel speed, U 00 = 2.24 m/s, the angle of attack was increased to 4° 
and the video camera was used to record the flow around the airfoil. A grabber computer 
card was used to extract the pictures from the recorded tape. 

A selected photo for the flow is shown in Fig. 4.3. They depict the development of 
the vortices around the airfoil and their shedding downstream into the wake. 

3. Biplane Configuration 

The biplane configuration was investigated for 3 different cases. The tunnel speed 
for all of these cases was U,^ = 2.24 m/s and Re=10,000. 
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Figure 4.3. Vortex Shedding Around NACA0014 at a = 4°, U 0 0 = 2.24 m/s, Re=10000 

a. Casel: minimum distance between the wings 

In this case, the two wings were set such that the distance between them was 
the minimum distance. The video camera was used to capture the presented consecutive 
frames with time step 1/30 of a second. These pictures are shown in Fig. 4.4, and Fig. 4.5. 

b. Case2: Maximum distance between the wings 

In this investigation, the distance between the two flapping wings was set to 
the maximum value. The wind tunnel speed was kept the same as for the previous case. A 
selected photo is presented in Fig. 4.6 

c. Cased: Wing flapping at k=l, f=5.6 Hz 

The model was set to flap at f=5.6 Hz, corresponding to a reduced frequency, 
k=l. A series of consecutive frames is presented to show the flow behavior during two 
flapping cycles. The first cycle is presented in Fig. 4.7, Fig. 4.8, and Fig. 4.9. The second 
cycle is presented in Fig. 4.10 and Fig. 4.11. 
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Figure 4.4. (a) Vortex Shedding Around Biplane with Minimum Distance between the Two 
Wings, Re=10000. The Time Step Between Frames =1/30 sec 
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Figure 4.5. (b) Vortex Shedding Around Biplane with Minimum Distance between the Two 
Wings, Re=10000, with time step=l/30 sec 
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Figure 4.6. Vortex Shedding Around Biplane with Maximum Distance between the Wings, 
Re=10000 

C. RESULTS 

During the experiment the strobe light along with the frame speed for the video 
camera was altered to be able to capture the flow details around the airfoil. When the 
frame speed was set to 1/60 of a second , it was impossible to capture any details, only 
the average stream lines were captured. When the frame speed was increased to 1/500 
of a second, the camera was able to capture the flow details. Hence, the strobe light was 
used to visually freeze the flow. As a result, the vortices that are shed downstream start to 
appear stationary, the strobe light frequencies were recorded and the shedding frequency 
was calculated as shown in Table 4.1 


52 




Figure 4.7. Vortex Shedding Around Biplane During Flapping, = 0° and $ = 30°, k=l, 
U=2.24 m/s, Re=10000 
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Figure 4.8. Vortex Shedding Around Biplane During Flapping, $ = 135° and $ = 180°, 
k=l, U=2.24 m/s, Re=10000 
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Figure 4.9. Vortex Shedding Around Biplane During Flapping, $ = 225° and = 360°, 
k=l, U=2.24 m/s, Re=10000 
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Figure 4.10. Vortex Shedding Around Biplane During Flapping, 4> = 0°, $ = 45° and 
4> = 135°, k=l, U=2.24 m/s, Re=10000 
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Figure 4.11. Vortex Shedding Around Biplane During Flapping, $ = 180°, $ = 225° and 
4> = 270°, k=l, U=2.24 m/s, Re=10000 
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Table 4.1. Shedding Frequency for Flapping Wing in Biplane Configuration 


Reynolds Number 

a [deg] 

Shedding Frequency, Hz 

5000 

0 

10.4 

5000 

4 

9.7 

7500 

0 

12.8 

7500 

4 

8.14 

10000 

0 

14.66 

10000 

4 

8.0 

15000 

0 

13.49 
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V. EXPERIMENTAL RESULTS 


A. THRUST MEASUREMENTS 

Jones [13] already provided extensive thrust measurements for the biplane configu¬ 
ration and several micro-air vehicle configurations. Therefore, in the present work, thrust 
generation by single plunging wings only was investigated. The effect of changing the 
flapping frequency on the thrust value was investigated. The effect of changing Reynolds 
number on the thrust force was investigated as well. The model was set to flap at frequen¬ 
cies between 2 to 7[Hz] with increment of l[Hz]. The wind tunnel speed was varied from 3 
to 9.5[m/s] with increment of 0.5[m/s]. The maximum flapping amplitude was 0.4C. 

The effect of increasing the flapping frequency on the thrust coefficient is illustrated 
in Fig. 5.1 for different Reynolds numbers. The error involved in the thrust measurements 
was included as error bars. 

In Fig. 5.2, the effect of varying Reynolds number on thrust coefficient is illustrated. 
The error in thrust measurement is included as error bars. 

B. STEADY FLOW MEASUREMENTS 

It was important to know detailed information about the flow behavior around the 
wing at the low speed of interest. A wing with a NACA0014 airfoil section was attached to a 
mechanism which allows a variable angle of attack. The LDV system was used to investigate 
the flow around the airfoil for three different cases. At zero angle of attack, the flow was 
probed at a free stream speed of 2 and 5 m/s, while for the third case the wing was at 4° 
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Figure 5.1. Thrust coefficient for single flapping wing in pure plunge with error included as 
error bar, h=0.4C, 
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Figure 5.2. Thrust coefficient versus reduced frequency, k, for single flapping wing in pure 
plunge, with measurements error included as error bar, h=0.4C 
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angle of attack and the free stream velocity was 2 m/s. For each case, 10 axial locations 
were investigated to cover from quarter chord upstream to two chords downstream of the 
trailing edge. At each axial location, 11 locations were probed in the vertical direction. 

1. First Case: U=2 m/s, a = 0° 

The measured axial velocity is non-dimensionalized and is illustrated as a contour 
and carpet plot in Fig. 5.3. The velocity profiles at each axial location are presented in 
Fig. 5.4 while the velocity profile over the entire domain is shown in Fig. 5.5. The measured 
turbulence intensity in the free stream was 2.8906. The turbulence intensity was normalized 
by this value and presented in Fig. 5.6 as carpet and contour plots. 

2. Second Case: U=5 m/s, a = 0° 

The non-dimensional axial velocity component u/U is shown in Fig. 5.7 as a contour 
plot and as a velocity profile. Also, the non-dimensional axial velocity component u/U at 
different axial locations X/C is presented in Fig. 5.8. The turbulence intensity at the free 
stream measured and found to be 2.6929. This value was used to normalize the measured 
turbulence intensities and presented in Fig. 5.9. 

3. Third Case: U=2 m/s, a = 4° 

For this case, the angle of attack was increased to 4 degrees. The flow was in¬ 
vestigated at 9 axial locations. At each axial location, 27 vertical locations were probed 
such that, —1 < Z/C < 1. The velocity values are shown as a carpet and contour plot in 
Fig. 5.10. The axial velocity profiles are presented in Fig. 5.11. The velocity profiles at 
different axial location are given in Fig. 5.12. 
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Figure 5.3. Axial Velocity at U=2 m/s and a = 0° 
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Figure 5.4. Velocity Profile at Different Axial Locations over NACA0014, U=2 m/s and 


a = 0 ° 
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Figure 5.5. Velocity Profile over NACA0014 at U=2 m/s and a = 0° 
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Figure 5.6. Nondimensional Turbulence Intensity over NACA0014 at U=2 m/s and a = 0° 
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Figure 5.7. Nondimensional Axial Velocity Component at U=5 m/s and a = 0° 
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Figure 5.8. Nondimensional Axial Velocity Component at U=5 m/s and a = 0° 
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Figure 5.9. Non dimensional Turbulence Intensity at U=5 m/s and a = 0° 
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The turbulence intensity was normalized by the measured turbulence intensity at 


the free stream which was 2.6929 and is presented in Fig. 5.13. 

C. UNSTEADY FLOW MEASUREMENTS 

The flapping model was used as the test model for the unsteady measurements. 
LDV was used along with the RMR to measure the axial velocity component and stamp it 
with the wing position through the triggering signal generated by the micro-switch. 

1. Biplane Configuration 

The flow around the biplane configuration during flapping was investigated. The 
wind tunnel speed was 2 m/s and the flapping frequency was 5 Hz, which corresponds 
to a reduced frequency, k, of 1. The Reynolds number was of 8763. Nine locations were 
investigated along the airfoil in the downstream direction. At each location, 17 points were 
measured. In total, 153 points were measured. At each point, the LDV was set up to obtain 
10000 samples or 8 minutes measuring time, whichever happens faster. It was found that 
10000 samples were enough to give acceptable statistics. 

The area investigated was 3.25 chord lengths in the downstream direction and 4.2 
chord lengths in the vertical direction. The locations investigated were X/C=-0.25, 0, 0.25, 
0.5, 0.75, 1.0, 1.5, 2.0 and 3.0. At each location, 17 points were investigated such that point 
number 9 was located nrid-way between the two flapping wing. 

The measured axial velocity component, u, at each position Z/C is presented as 
a function of the plunge phase angle, $. In the same time, the measured u is presented 
as carpet and contour plots at each axial location. Nondimensional turbulence intensity, 
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Figure 5.10. Axial Velocity at U 00 = 2 m/s and a = 4° 
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Figure 5.11. Axial Velocity Profile at Measured Axial Locations, U 00 = 2m/s and a = 4° 
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Figure 5.12. Axial Velocity Profile at Uoo = 2 m/s and a = 4° 
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Figure 5.13. Normalized Turbulence Intensity at Uoo = 2m/s and a = 4° 
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Table 5.1. Biplane Configuration 


X/C 

u/U ooFs.T 

u/U as carpet and contour plots 

Ti/Uoo 

-0.25 

Fig. B.l 

Fig. 5.14 

Fig. 5.15 

0.0 

Fig. B.2 

Fig. 5.16 

Fig. 5.17 

0.25 

Fig. B.3 

Fig. 5.18 

Fig. 5.19 

0.5 

Fig. B.4 

Fig. 5.20 

Fig. 5.21 

0.75 

Fig. B.5 

Fig. 5.22 

Fig. 5.23 

1.0 

Fig. B.6 

Fig. 5.24 

Fig. 5.25 

1.5 

Fig. B.7 

Fig. 5.26 

Fig. 5.27 

2.0 

Fig. B.8 

Fig. 5.28 

Fig. 5.29 

3.0 

Fig. B.9 

Fig. 5.30 



Ti/U 0 o is presented as carpet plot at each location. The related figures for each axial 
position X/C is presented in Table 5.1 

The two wings start flapping above the middle point for each one of them at plunge 
phase angle, $ = 0 in opposite directions. This means that the upstroke starts at $ = 0. 

a. Average Velocity 

The data extracted from the LDV measurements was used to calculate the 
average velocity profile at each location. The obtained velocity profiles are given in Fig. 5.32. 
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Biplane configuration, X/C = - 0.25, k=1, Re=8760 
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Figure 5.14. Measured Axial Velocity Component, u, at X/C=-0.25 


76 





























































Figure 5.15. Nondimensional turbulence intensity at X/C=-0.25 
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Biplane configuration, X/C=0, k=1, Re=8760 



Biplane configuration, X/C=0, k=1, Re=8760 



Figure 5.16. Measured Axial Velocity Component, u, at X/C=0, the Leading Edge 
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Figure 5.17. Nondimensional Turbulence Intensity at the Leading Edge 
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Biplane configuration, X/C=0.25, k=1, Re=8760 
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Figure 5.18. Measured Axial Velocity Component, u, at X/C=0.25 
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Figure 5.19. Nondimensional Turbulence Intensity at X/C=0.25 
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Biplane configuration, X/C=0.5, k=1, Re=8760 
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Figure 5.20. The Measured Velocity over the Complete Domain at X/C = 0.5 
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Nondimensional Turbulence Intensity 
Biplane configuration, X/C=0.5, k=1, Re=8760 
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Figure 5.21. Nondimensional Turbulence Intensity, Ti/Uoo at X/C = 0.5 
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Biplane configuration, X/C=0.75, k=1, Re=8760 
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Figure 5.22. The Measured Axial Velocity Component over the Complete Domain at 


X/C=0.75 
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Nondimensional Turbulence Intensity 
Biplane configuration, X/C=0.75, k=1, Re=8760 
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Nondimensional Turbulence Intensity 
Biplane configuration, X/C=0.75, k=1, Re=8760 
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Figure 5.23. Nondimensional Turbulence Intensity at X/C=0.75 
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Biplane configuration, X/C=1, k=1, Re=8760 

u/U M 

■ 1.72 
= 1.51 
= 1.29 
= 1.07 

- 0.85 
I 0.64 

■ 0.42 

3 0.20 


0 90 180 270 360 

Plunge phase angle, [deg] 


Biplane configuration, X/C=1, k=1, Re=8760 




Figure 5.24. Measured Velocity over the Complete Domain 
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Nondimensional Turbulence Intensity 
Biplane configuration, X/C=1, k=1, Re=8760 
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Figure 5.25. Nondimensional Turbulence Intensity 
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Biplane configuration, X/C=1.5, k=1, Re=8760 
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Figure 5.26. Measured Velocity over the Complete Domain at X/C = 1.5 


88 



































Nondimensional Turbulence Intensity 
Biplane configuration, X/C=1.5, k=1, Re=8760 
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Figure 5.27. Nondimensional Turbulence Intensity at X/C = 1.5 
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Biplane configuration, X/C = 2, K=1, Re=8760 
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Figure 5.28. Measured Velocity over the Complete Domain at X/C=2.0 
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Figure 5.29. Nondimensional Turbulence Intensity 
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Biplane configuration, X/C=3, k=1, Re=8760 
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Figure 5.30. Measured Velocity over the Complete Domain at X/C = 3.0 
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Figure 5.33. Flow Around NACA0014 in the Upper Flapping Stroke for Biplane Configu¬ 
ration 


b. Summary 


The measured instantaneous velocities were used to obtain the flow field 


over the whole airfoil at different angles, 4>, which helps to understand the flapping motion. 


The flow around the airfoil is illustrated over a complete flapping cycle in Fig. 5.33, and 


Fig. 5.34. 
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Figure 5.34. Flow Around NACA0014 in the Lower Flapping Stroke for Biplane Configu¬ 
ration 
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2 . 


Wing in Ground Effect 


For this experiment, the flapping model was configured to use only one wing. A brass 
block with the same weight as the removed wing was added to the model to ensure stability. 
A fixed flat plate was mounted underneath the wing such that the distance between the 
plate and the flapping wing in the lowest point in the downstroke was 22 mm. 

At first, the boundary layer over the flat plate was studied. The velocity profile was 
measured at two different locations,namely at X/C = 5 and X/C = 13, where X is the 
distance measured from the leading edge of the flat plate, and C is the airfoil chord. The 
velocity profile at X/C = 5 was measured twice, once with the flapping wing at the lowest 
point of the downstroke , while the other was at the highest point of the upstroke. For 
X/C = 13, the velocity profile was measured only once when the wing was at the highest 
point. 

The velocity profile is shown in Fig. 5.35. The standard deviation is presented as 
error bars. The turbulence intensity is presented in the same figure. The distance between 
the ground plane and the flapping-wing mid-point was Z/C = 0.84. This value was selected 
such that the wing will be immersed in the boundary layer in the downstroke part of the 
flapping cycle. 

After the boundary layer was investigated, the model was set to flap at 5 Hz which 
corresponds to a reduced frequency, k = 1, at the used tunnel speed of 2 m/s. The 
Reynolds number therefore was 8874. The flow around the flapping airfoil was probed in 80 
points by LDV, distributed as a set of ten points at eight different axial locations, namely, 
X/C = —0.25,0,0.25,0.5,0.75,1.0,1.5,2.0. At each location, the flow was investigated in 
the vertical direction between —0.69 < Z/C < 1.4. The closest point of measurements to 
the ground plane was at a distance of 3 mm. 
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Figure 5.35. Boundary Layer over the Ground Plane 
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Table 5.2. Flapping Wing in Ground Effect 


X/C 

u/UooVs .$ 

u/U as carpet plot 

-0.25 

Fig. D.l 

Fig. 5.36 

0 

Fig. D.2 

Fig. 5.37 

0.25 

Fig. D.3 

Fig. 5.38 

0.50 

Fig. D.4 

Fig. 5.39 

0.75 

Fig. D.5 

Fig. 5.40 

1.0 

Fig. D.6 

Fig. 5.41 

1.5 

Fig. D.7 

Fig. 5.42 

2.0 

Fig. D.8 

Fig. 5.43 


The axial velocity component, u, was measured as a function of the encoder position 
which is related to the plunge phase angle at each probed point. The results are presented 
in the subsequent figures shown in table 5.2 

The obtained measurements was integrated together to build a three dimensional 
block of the nondimensional velocity u/U inf. The integration of the measurements in 3D 
block with X/C as the first coordinates, Z/C as the second coordinates and Plunge Phase 
Angle as the third coordinates, allows information extraction in different useful directions. 
Eight slices were taken at different plunge phase angle, starting at zero and by incremental 
step of 45 degree to illustrate the changes in the velocity distribution around the airfoil in 
both upstroke and downstroke, which is illustrated in Fig. 5.44 and Fig. 5.45. 
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Figure 5.36. Nondimensional Axial Velocity at X/C=-0.25 for Flapping Wing in Presence 


of Ground Effect 


98 












































































































Flapping wing in presence of ground plane 
X/C=0, k=1, Re=8760 
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Figure 5.37. Flapping Wing in Presence of Ground Effect, X/C=0 


99 































































Flapping wing in presence of ground plane 
X/C=0.25, k=1, Re=8760 
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Figure 5.38. Flapping Wing in Presence of Ground Effect, X/C=0.25 
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Flapping wing in presence of ground plane 
X/C=0.5, k=1, Re=8760 
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Figure 5.39. Flapping Wing in Presence of Ground Effect, X/C=0.5 
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Flapping wing in presence of ground plane 
X/C=0.75, k=1, Re=8760 
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Figure 5.40. Flapping Wing in Presence of Ground Effect, X/C=0.75 
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Flapping wing in presence of ground plane 
X/C=1, k=1, Re=8760 
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Figure 5.41. Flapping Wing in Presence of Ground Effect, X/C=l. 
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Flapping wing in presence of ground plane 
X/C = 1.5, k=1, Re=8760 
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Figure 5.42. Flapping Wing in Presence of Ground Effect, X/C=1.5. 
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Flapping wing in presence of ground plane 
X/C=2, k=1, Re=8760 
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Figure 5.43. Flapping Wing in Presence of Ground Effect at X/C=2.0 
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Figure 5.44. Flapping Wing in Ground Plane During Upstroke. 
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Figure 5.45. Flapping Wing in Ground Plane During Downstroke. 
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Table 5.3. Single Flapping Wing 


X/C 

u/UooVs .<& 

u/U as carpet plot 

-0.2 

Fig. C.l 

Fig. 5.49 

0 

Fig. C.2 

Fig. 5.50 

0.2 

Fig. C.3 

Fig. 5.51 

0.6 

Fig. C.4 

Fig. 5.52 

1.0 

Fig. C.5 

Fig. 5.53 


3. Single Wing Configuration 

In the third experiment a single sinusoidally plunging wing was investigated at a 
reduced frequency of k = 1.0. The same wing previously used for the ground effect was 
set to flap at a frequency of 5Hz, while the wind tunnel speed was set to 2 m/s. The LDV 
along with the RMR was used to measure the velocity at 5 locations in the axial direction. 
Namely at X/C=-0.2, 0, 0.2, 0.6 and 1.0. At each location, 15 points were investigated in the 
vertical direction to cover —1.4 < Z/C < 1.4. The flapping range was —0.4 < Z/C < 0.4. 

The non-dimensional axial velocity distribution is shown in Fig. 5.46 and Fig. 5.47. 
At each position X/C , the instantaneous velocity is given for each measurement point Z/C 
as mentioned in table 5.3. 

The average velocity profile around the airfoil is shown in Fig. 5.48. It is presented 
at the axial locations X/C to show how the variation occurs in the flow around the airfoil. 
The non-dimensional turbulence intensity around the airfoil is represented as well. The 
turbulence intensity at the free stream was used to non-dimensionalize the measured values. 


108 





Figure 5.46. Single Flapping Wing During Upstroke, k=l, u=2[m/s], f=5[Hz], Re=8760 
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Figure 5.47. Single Flapping Wing During Downstroke, k=l, u=2[m/s], f=5[Hz], Re=8760 
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Figure 5.48. Average Velocity Profile and the Associated Turbulence Intensity Around 

Single Flapping Wing, k=l, u=2 m/s, f=5 Hz 
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Figure 5.49. Single Flapping Wing at X/C=-0.2, k=l, u=2[m/s], f=5[Hz], Re=8760 
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Pure plunge, k=1, Re=8760, X/C = 0 
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Figure 5.50. Single Flapping Wing at X/C=0.0,k=l, u=2[m/s], f=5[Hz], Re=8760 
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Figure 5.51. Single Flapping Wing at X/C=0.2, k=l, u=2[m/s], f=5[Hz], Re=8760 
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Pure Plunge, k=1, Re=8760, X/C = 0.6 
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Figure 5.52. Single Flapping Wing at X/C = 0.6, k=l, u=2 m/s, f=5 Hz, Re=8760 
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Figure 5.53. Single Flapping Wing at X/C = 1.0, k=l, u=2m/s, f=5 Hz, Re=8760 
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VI. ANALYSIS OF THE EXPERIMENTAL 


RESULTS 


A. STEADY FLOW MEASUREMENTS 


In the steady flow investigation, the axial velocity component, u, and the turbulence 
intensity were measured at 2 m/s freestream velocity at two different angles of attack, 
namely 0° and 4°. It was also measured at 5 m/s for zero degree angle of attack. 

The obtained results are presented in Fig. 5.3 through Fig. 5.12. For the three cases 
the turbulence intensity has the highest value in the wake, as expected. Also, a high velocity 
area can be seen on the upper and lower surfaces for the airfoil at zero angle of attack for 
U=2 and 5 m/s, Fig. 5.3 and Fig. ??. Meanwhile, for the case of angle of attack =4° the 
highest value for the velocity can be seen on the upper surface only, as shown in Fig. 5.10. 

The normalized velocity profiles, u/Uoo are presented in Fig. 5.4, Fig. 5.8 and 
Fig. 5.12. For the case of 4 degrees angle of attack, the changes in the wake behind the 
airfoil are clearly distinguished. For the case of Uoo = 2 and AOA=0°, the highest value of 
the turbulence intensity is found on the lower surface near the mid chord position. Also, 
two areas with high velocity and high turbulence intensity can be seen at one chord behind 
the trailing edge in Fig. 5.6. 
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B. THRUST MEASUREMENTS 


The thrust generation by wing flapping was investigated as well. The thrust co¬ 
efficient as function of flapping frequency and reduced frequency is presented for different 
Reynolds numbers in Fig. 5.1 through Fig. ??. It can be seen that as the reduced frequency 
increases, the thrust coefficient increases for the same Reynolds number. Also, for the same 
reduced frequency, the thrust coefficient increases as the Reynolds number increases. 

C. FLAPPING BIPLANE FLOW MEASUREMENTS 

The flow between and around the two flapping airfoils was investigated. The axial 
velocity component, u, is presented over the whole domain as contour and carpet plots to 
present the changes during a complete flapping cycle at each location X/C. The same for 
the turbulence intensity. Also, at each location, X/C , the measured velocity component u 
is presented as nondimensional value at different vertical locations Z/C. It is important to 
note that in some plots, like in Fig. B.2, some velocity measurements yield zero velocity 
values. This is due to the fact that the LDV system was unable to detect any scattered 
light from the particles passing through the laser probe. This might be due to the fact that 
the measurements are made in high vortical regions with too few seed particles. It might 
also be due to laser beam blockage as the wings pass through the beams during the flapping 
cycle. 

All the velocity measurements were processed and the velocity distributions around 
the airfoil at different plunge phase angles were extracted and presented in Fig. ??, and 
Fig. ?? for a complete cycle. Since the flow around the two airfoils is antisymmetric around 
the plane of symmetry, the flow around the upper airfoil will only be described here. Further- 
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more, it is important to recognize the difficulty in identifying the precise flow characteristics 


based on the axial flow field information only. Therefore, in the following paragraphs we —- 
describe the locations of high and low axial velocities —- of distinct vortices (which could 
be better identified with u- and v-velocity information). 


At plunge phase angle, $ = 0°, the upper airfoil is located at ZjC = 1.1. It is 
noticed that the velocity near the upper surface is higher than that near the lower surface. 
The wake area behind the airfoil has a high velocity value. As the wing moves down to 
$ = 45°, the velocity increases on the upper surface while the area of high velocity in the 
wake is moving downstream. As the wing reaches the mean position at $ = 90°, ZfC = 0.7, 
the high velocity region still exists on the upper surface but moves downstream towards the 
trailing edge. An area of low velocity is noticed on the lower surface. The area of high 
velocity in the wake is moving further downstream. As the wing keeps moving down during 
the downstroke, $ = 135°, ZfC = 0.5, the high velocity region on the upper surface is 
moving further downstream. The velocity underneath the wing is lower than that above it. 
A low velocity area is noticed in the wake behind and between the two airfoils. Meanwhile, 
the velocity between the two airfoils starts to increase, especially within the first quarter 
chord. The increase in the axial velocity along the centerline between the two airfoils can 
also be seen quite clearly in Fig. 5.16, Fig. 5.18 and Fig. 5.20. The highest velocities occur 
when the two airfoils are close together around $ = 180°. The wing reaches the lowest 
point during the downstroke at, $ = 180°, Z/C 1 = 0.3. At this location, it appears that 
three vortices exist on the upper surface. For two of them, the axial velocity component 
is small while the third one has a higher value. Comparing with the computational results 
and the flow visualization, we can conclude that there are two spots with low axial velocity 
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because we have two clockwise vortices which decreases the axial velocity. The opposite 


situation happens with the third vortex which is a counterclockwise vortex. 


As the wing starts the upstroke, $ = 225°, Z/C = 0.5, the effective angle of attack 
changes direction and becomes negative. Hence, the velocity starts to increase and build up 
on the lower surface while the upper surface has a lower velocity. As the airfoil reaches the 
mean position during the upstroke, $ = 270°, Z/C = 0.7, the high velocity area between 
the two airfoils starts to move downstream to almost one chord behind the trailing edge. 
As the airfoil keeps moving up during the upstroke at $ = 315°, Z/C = 0.9, the velocity at 
the lower surface is higher than that on the upper surface. Three vortices can be noticed. 
Two counterclockwise and one clockwise. The high velocity area between the two airfoils 
now is moving downstream. 


The axial velocity component, u, was measured over a complete cycle at several 
positions X/C. These instantaneous velocity values were averaged to obtain a velocity 
profile at each axial location, X/C and presented in Fig. ??. At the leading edge, it is 
seen that the velocity decreases dramatically within the flapping margin in such a way that 
there are two peaks in the mean flapping position of the two airfoils. Meanwhile, between 
the leading and trailing edges the velocity is higher underneath the upper flapping margin 
than the lower part of it. In the wake behind the airfoil, two increases in velocity can be 
noticed behind the trailing edge which indicate the shed vortices. 
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D. FLOW MEASUREMENTS FOR FLAPPING AIRFOIL 


IN GROUND EFFECT 

As the wing starts flapping from the uppermost position it is noticed that the velocity 
starts to increase between the wing and the ground plane. As the wing keeps moving down 
during the downstroke to = 45°, it is noticed that an area of high velocity develops near 
the lower surface of the airfoil downstream of the midchord position. As the wing moves 
toward the mean position, $ = 90°, two distinct areas can be seen between the wing and the 
ground plane. The first is a low velocity region just behind the leading edge. The second is 
a high velocity region near the trailing edge. 

As the wing keeps moving down to, $ = 135°, the high and low velocity regions 
between the wing and the ground plane move downstream. As the wing keeps moving 
down to T = 180°, a small area of higher velocity develops on the upper surface close to 
the leading edge. Meanwhile, the low velocity region keeps moving downstream near the 
ground plane. 

During the upstroke at = 225°, the effective angle of attack changes direction, 
and the velocity near the upper surface becomes less than the free stream velocity. The 
velocity starts to increase a little bit on the lower surface. At = 270°, four regions can be 
distinguished , two with high velocity and two with low velocity. As the wing approaches 
= 315°, a low velocity region is noticed near the upper surface of the airfoil. 

E. SINGLE FLAPPING AIRFOIL FLOW MEASUREMENTS 

The wing starts the flapping cycle at $ = 0° from the uppermost position, Z/C = 

0.4. It is noticed that there is a high velocity region near the lower surface of the airfoil. 
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This vortex was generated during the previous upstroke. As the airfoil moves down, the 
effective angle of attack increases and the velocity increases on the upper surface. As the 
airfoil approaches the mean position, $ = 90°, the velocity near the upper surface becomes 
higher than the freestream. At the same time the velocity on the lower surface decreases 
and becomes less than the freestream value. 

At = 135°, the velocity on the upper surface is higher than the free stream, 
indicative of a vortex. At $ = 180°, the downstroke ends. At that position, a high velocity 
region is noticed on the upper surface of the airfoil. As the wing starts the upstroke, 
$ = 225°, the effective angle of attack changes sign and the velocity near the lower surface 
starts to increase and becomes higher than that on the upper surface. As the wing keeps 
moving upward, the velocity increases on the lower surface and the high velocity region 
starts to sweep backward until it is almost shed from the lower surface at $ = 315°. 
The comparison between the biplane measurements and the single airfoil results shows the 
following significant differences. 

A direct comparison between the biplane, flapping wing in presence of ground plane 
and single flapping wing at different axial locations X/C for different phase angles is pre¬ 
sented in Tables 6.1, 6.2, 6.3, 6.4, 6.5 and 6.6. The axial velocities listed in these 
tables were measured at Z/C=-0.7 from the mean flapping position, which is the plane of 
symmetry for the biplane case. It is readily seen that the biplane configuration, in general, 
generates larger disturbances than the single airfoil or the airfoil in ground effect, thus giving 
further quantitative evidence for the greater thrust generation capability of the biplane. 
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Table 6.1. Measured Axial Velocity Component, u/Uoo, at Different Plunge Phase Angle, 
$, Z/C = -0.7, k=l, Re=8760, h=0.4, X/C=-0.25 


Plunge phase angle 

$ [deg] 

Biplane 

X/C=-0.25 

Ground plane 

X/C=-0.25 

Single flapping wing 

X/C=-0.2 

0.83 

0.863769619 

0.750694541 

0.9834992 

45.83 

0.698112989 

0.608992773 

0.936257251 

90.83 

0.804731722 

0.741699892 

0.961867373 

135.83 

1.087860634 

0.962276964 

1.034629426 

180.83 

1.220980759 

1.067223162 

1.075115023 

225.83 

1.244995684 

1.08709762 

1.073714743 

270.83 

1.175133207 

0.995445301 

1.055003501 

315.83 

1.039729395 

0.909875855 

1.018581216 
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Table 6.2. Measured Axial Velocity Component, u/Uoo, at Different Plunge Phase Angle, 
$, Z/C = —0.7, k=l, Re=8760, h=0.4, X/C=0, Leading Edge 


Plunge phase angle 

<f> [deg] 

Biplane 

x/c=o 

Ground plane 

x/c=o 

Single flapping wing 

x/c=o 

0.83 

0.811684173 

0.629574881 

0.972944589 

45.83 

0.569932001 

0.433633243 

0.939377876 

90.83 

0 

0.33616 

0.966178236 

135.83 

1.13911082 

1.022446754 

1.033916783 

180.83 

1.261565692 

1.114126203 

1.077715543 

225.83 

1.23081558 

1.036677854 

1.085364573 

270.83 

1.164183097 

0.976808788 

1.065453091 

315.83 

1.020320469 

0.883181446 

1.027738048 
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Table 6.3. Measured Axial Velocity Component, u/Uoo , at Different Plunge Phase Angle, 
$, Z/C = -0.7, k=l, Re=8760, h=0.4, X/C=0.25 


Plunge phase angle 

<f> [deg] 

Biplane 

X/C=0.25 

Ground plane 

X/C=0.25 

Single flapping wing 

X/C=0.2 

0.83 

0.774816871 

0.802111299 

1.004863473 

45.83 

0.575344291 

1.245323812 

0.946941888 

90.83 

0 

0.312615508 

0.985342068 

135.83 

1.256438259 

0.861304117 

1.066218244 

180.83 

1.301228459 

1.016753033 

1.105781156 

225.83 

1.255636796 

1.062136096 

1.110977195 

270.83 

1.179536426 

1.050574114 

1.098187137 

315.83 

1.024004303 

0.970808273 

1.036912382 
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Table 6.4. Measured Axial Velocity Component, u/Uoo, at Different Plunge Phase Angle, 
$, Z/C = -0.7, k=l, Re=8760, h=0.4, X/C=0.5 


Plunge phase angle 

$ [deg] 

Biplane 

X/C=0.5 

Ground plane 

X/C=0.5 

Single flapping wing 

X/C=0.6 

0.83 

0.900424486 

0.904803634 

1.01315013 

45.83 

0.734425061 

1.263628183 

0.95585 

90.83 

0 

0.97328 

0.937267453 

135.83 

1.182476733 

0.978371007 

0.98485 

180.83 

1.269807243 

1.068632203 

1.0282 

225.83 

1.275504391 

1.091999802 

1.050573825 

270.83 

1.213844841 

1.079994272 

1.05694942 

315.83 

1.093128082 

1.014842246 

1.038895279 
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Table 6.5. Measured Axial Velocity Component, u/Uoo, at Different Plunge Phase Angle, 
$, Z/C = -0.7, k=l, Re=8760, h=0.4, X/C=0.75 


Plunge phase angle 

$ [deg] 

Biplane 

X/C=0.75 

Ground plane 

X/c=0.75 

0.83 

0.981739194 

0.937109564 

45.83 

0.97140708 

0.805147089 

90.83 

0.438033383 

2.1488 

135.83 

0.905909197 

0.715965832 

180.83 

1.167707603 

0.891507233 

225.83 

1.273553842 

1.016394608 

270.83 

1.230279662 

0.993620106 

315.83 

1.209224358 

1.010877598 
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Table 6.6. Measured Axial Velocity Component, u/Uoo, at Different Plunge Phase Angle, 
$, Z/C = —0.7, k=l, Re=8760, h=0.4, X/C=l, Trailing Edge 


Plunge phase angle 

<f> [deg] 

Biplane 

X/C=l 

Ground plane 

X/C=l 

Single flapping wing 

X/C=l 

0.83 

1.132776365 

1.008925295 

1.037340368 

45.83 

1.068528958 

0.902613539 

1.023719744 

90.83 

0.485672157 

1.128697907 

0.943016103 

135.83 

0.665957903 

0.49128 

0.936744249 

180.83 

1.073062535 

0.763821508 

0.96765103 

225.83 

1.216060934 

0.935216116 

1.027968094 

270.83 

1.286237237 

1.015900262 

1.049592418 

315.83 

1.267735992 

1.028645633 

1.066170734 
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VII. CONCLUSIONS AND 


RECOMMENDATIONS 

A wind tunnel investigation of low speed flow over flapping airfoils and airfoil com¬ 
binations was performed using flow visualization and laser doppler velocimetry. Specifically, 
three cases were studied: A NACA0014 airfoil oscillating in a sinusoidal plunge mode, A 
NACA0014 airfoil oscillating in a sinusoidal plunge mode near a ground plane, and two 
NACA0014 airfoils arranged in a biplane configuration and oscillating in counterphase in a 
sinusoidal plunge mode. The plunge amplitude-to-airfoil chord ratio was 0.4, the reduced 
frequency of oscillation was 1.0 and the Reynolds number based on airfoil chord was set at 
8760. 

Conditionally sampled measurements of the axial flow velocity were taken at nu¬ 
merous flow field points providing detailed information about the flow features generated 
by this type of flapping motion. These measurements were complemented by tinre-averaged 
flow field data and by visualization of the instantaneous flow field at various points during 
the flapping cycle. Furthermore, the thrust generated by the sinusoidal plunge motion was 
measured with a laser range finder. The results shows that vortex shedding occurs both 
from the airfoil leading and trailing edge. 

This experimental evidence for the occurrence of dynamic stall at the chosen test 
conditions of h=0.4, k=l and Reynolds nunrber=8760 provides credence to the Navier- 
Stokes predictions summarized in Chapter II. Dynamic stall may well be an important, and 
even critical, ingredient to achieve optimum thrust and lift on low Reynolds number flight 
vehicles, such as birds, and micro-air vehicles. 
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However, before more specific conclusions can be stated about the role of dynamic 


stall for effective thrust and lift generation, the following additional studies should be pur¬ 
sued: 

1. Repeat the measurements reported in this thesis with a Particle Image Velocimetery 
PIV system so that complete flow field information becomes available. 

2. Extend the range of test parameters so that the dynamic stall boundary for sinu¬ 
soidal plunge oscillation can be determined. This requires a repetition of the tests 
for different values of plunge amplitude and frequency and tunnel speeds. Also, it 
is recommended to test for plunge amplitudes greater than 0.4 so that the change 
in flow features can be measured as the airfoil is oscillating in deep dynamic stall. 
Unfortunately, this requires a significant modification of the model which is currently 
limited to h=0.4 

3. Extend the tests to explore the effect of combined pitch and plunge motions. As 
shown by the potential flow predictions of Jones and Platzer [13] and the Navier- 
Stokes predictions of Isogai et al [11] the phase angle between pitch and plunge and 
the amplitude ratio between pitch and plunge are two additional critical parameters. 
Detailed experimental flow investigation is needed to verify these predictions. 

4. Detailed comparisons between the Navier-Stokes predictions of Castro [5] and the 
present experiments for sinusoidally plunging airfoils need to be made so that the 
predictive capability of the Navier-Stokes computations can be assessed. 

5. If the Navier-Stokes predictions are found to be reasonably accurate, additional com¬ 
putations are recommended to identify the parameter combination for optimum thrust 
and lift generation 
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APPENDIX A: CALIBRATION OF WIND 


TUNNEL 


At the beginning of the experimental work, the wind tunnel test section was cal¬ 
ibrated to document the flow quality. The calibration was done in three steps. The first 
step was a survey of the turbulence characteristics over the whole wind tunnel speed range. 
In the second step, a pitot tube was used for a velocity survey. The third step consisted of 
a detailed survey using LDV in one vertical plane located one half chord length upstream 
of the leading edge location of the wing in all experiments. 


A. Turbulence Characteristics 


LDV was used to measure the axial velocity component, u, at a fixed point. It 
was intended to take such measurements at the exact center of the test section but due to 
traverse mechanism limitations, this point was located at 68 cm height from the floor and 
71 cm from the side wall. This is only 7 and 4 cm, respectively from the actual center of 
the test section. The laser probe was aimed to be located 43 cm underneath the pitot tube. 
The axial velocity component, u, turbulence intensity, standard deviation and the number 
of particles passing through the probe volume, were recorded and are presented in Table 
A.l. 
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Table A.l. Wind Tunnel Calibration 


U mean 

Standard deviation 

Turbulence Intensity 

Number of particles 

1.06132 

0.0403795 

3.80464 

2965 

2.00892 

0.0476697 

2.3729 

2983 

3.00066 

0.0808584 

2.69469 

2979 

3.96654 

0.103202 

2.60182 

2997 

5.12694 

0.134258 

2.61867 

2990 

6.02697 

0.158603 

2.63155 

2993 

7.04868 

0.193352 

2.7431 

2980 

7.96251 

0.215616 

2.70789 

2990 


B. Pitot Tube Survey 

In his thesis, Lund [21], refurbished the wind tunnel. He calibrated the wind tunnel 
using both a pitot tube and LDV and presented an empirical formula (A.l) which relates 
the pitot tube output voltage to the LDV measurements. He did two sets of calibration. In 
the first one, the LDV and pitot tube were co-located, with the LDV probe 0.1 [m] above 
and 0.05 [m] in front of the pitot tube. In the second experiment, the LDV probe was 2.4 
[m] downstream and 0.1 [m] lower. The empirical formula for the second case, (A.2), was 
close to the first one. 


V = 0.4819 x ( mVolt ) 0A89 


(A.l) 


V = 0.4593 x (mVolt) 0 ' 4979 


(A.2) 
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Where, V, is the velocity in m/s and mVolt is the output of the pitot tube system 
in millie Volts. 

In the present work, the pitot tube was moved 2.4 nr downstream from the location 
used by Lund. The new location was selected such that it is located above the test model. 
A few measurements were taken for the wind tunnel speed by LDV and compared with the 
detected pitot tube voltage. It turned out that equation (A.l) is still valid and so it was 
used in the following survey. 

In the present work, a Pitot static tube was used to investigate the flow velocity along 
the test section. The pitot tube was used to probe 168 locations distributed such that the 
test section could be surveyed completely. In the streamwise direction, X, 7 locations were 
investigated with 0.3 nr apart. At each X location, a vertical plane ( Y,Z ) was investigated. 
Along the lateral direction, Y, six points were investigated. The first and sixth point was 
located 0.15 nr from the side walls, while the distance between the points was 0.3 nr. In the 
vertical direction, Z, four locations were investigated, 0.3 nr apart and the first point was 
0.3 nr from the floor. The investigated points are shown in Fig. A.l. 

Since the point located at (0,0.76,0.61) was almost in the middle of the test section, 
the axial velocity obtained from this point was used as the axial component of the free 
stream velocity, U. Data obtained from the pitot tube are presented in Fig. A.2. At each 
X location 24 measurements obtained at the associated ( Y ., Z) plane are plotted. Obtained 
data at the same vertical distance Z have the same symbol. 

The main purpose of this experiment was to investigate the flow quality in the test 
section. The obtained velocity measurements are presented as a 3D contour graphs. The 
highest nondinrensional velocity recorded was 1.011 which was located at only one spot 
close to the side wall, far away from the wing location. It is seen that the non-dinrensional 
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Figure A.l. Locations Probed by The Pitot Tube Along The Wind Tunnel Test Section 
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Figure A.2. Nondimensional Axial Velocity Obtained by the Pitot Tube 
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velocity values range from 0.978 to 1.011. In Fig. A.3 the flow survey in a vertical planes 
perpendicular to the X axis is presented at different X locations. Also, a set of surveys 
taken perpendicular to Y and Z is shown in Fig. A.4 and Fig. A.5. 

C. LDV Survey 

The LDV system was used to scan two vertical planes. One (V, Z) and one ( Y ., Z ). 

For the (X,Z) flow survey, the axial velocity component, u, was measured at the 
streamwise stations X = —0.5 C, X = 0.5 C and X = 1.5 C. The measured axial velocity 
component, u, is presented versus the vertical distance, Z, in Fig. A.6. The turbulence 
intensity is illustrated in Fig. A.7. 

It is seen that the u velocity in the wind tunnel increases with increasing distance 
from the floor, effectively exposing the model to a slight shear flow. This becomes very 
obvious when the velocities are non-dimensionalized by the local free stream velocity which 
is taken at the center of the tunnel. Curve fitting produces the polynomial, Eq. (A.3), as 
analytical approximation for this shear flow, valid between 0.381 < Z < 0.508 [nr] 

u(z) = 11.069 * z 2 - 8.771 * z + 4.6758 (A.3) 

For the second flow survey, a (V, Z) plane located a half-chord upstream of the 
model’s leading edge was chosen. Nine points were used to scan the lateral direction Y 
from 0.508 nr to 0.7112 nr measured from the tunnel side wall which has the laser probe. 
The vertical direction Z was scanned with 18 points from 0.5472 nr to 0.6858 nr. In total, 
162 points were used to scan this plane. The scanned plane location for this flow survey is 
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Figure A.3. Nondimensional Axial Velocity Component Measured by Pitot Tube at Different 


Axial Locations 
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Figure A.4. Nondimensional Axial Velocity Component Measured by Pitot Tube at Different 


lateral Locations 
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Figure A.5. Nondimensional Axial Velocity Component Measured by Pitot Tube at Different 


Vertical Locations 
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Figure A.6. Measured Velocity 



Figure A.7. Turbulence Intensity % 
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shown in Fig. A.8. The X axis scale is inverted on purpose to make it easy for the reader 


to compare with the results obtained from the pitot tube measurements. 

The measured axial velocity, u, was non-dimensionalized by the measured value 
at Y = 0.7112 m and Z = 0.6858 m. This point was closest to the center of the test section 
which can be reached by the LDV probe. This location was used to measure the free stream 
velocity for all the experiments in the present work. 

The non-dimensional axial velocity u/Uoo is presented in Fig. A.9 as it varies with 
vertical distance. It is noticed that the velocity increases a little as one approaches the center 
of the tunnel. The wing was located at Z = 0.45m in this work. The standard deviation 
of the measurements was added as an error bar in the same figure. The non-dimensional 
axial velocity u/Uoo is shown also as a contour lines in Fig. A.10 The turbulence intensity 
of the measurements is presented in two ways, as scattered points and as contour lines, in 
Fig. A. 11 gives complete information about the whole domain. 

D. DATA SCATTER AND THE STANDARD DEVIA¬ 
TION FOR THE MEASUREMENTS 

Testing the validity of the measurements and its scatter was important. The axial 
velocity component, u, was measured by LDV at three different locations. The first location 
was at the far upstream almost at the middle of the test section. The second location was 
at one chord upstream of the leading edge. The third one was at one chord downstream of 
the trailing edge. The standard deviation for each set was calculated and shown in table 
A.2. Meanwhile, the measured velocity at the three locations is presented in Fig. A.12. 
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Figure A.8. The Scanned Plane 
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Figure A.9. u/Uoo 
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Figure A.10. u/Uoo 


Table A.2. Standard Deviation for the Measured Axial Velocity Component 


Location of measurements 

Standard deviation of measurements 

Far free stream 

0.01229 

One chord upstream of leading edge 

0.014044 

One chord downstream of the trailing edge 

0.010138 
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Figure A. 11. Turbulence Intensity 
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Figure A. 12. Measured Axial Velocity Component, u, at Different Locations. 
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APPENDIX B: UNSTEADY MEASUREMENTS 


FOR THE BIPLANE CASE 

1. First location, X/C = —0.25 

The measured velocity, u, is presented as function of the plunge phase angle, 4>, at 
all vertical locations Z/C in Fig. B.l. 

2. Second Location, X/C = 0, (leading edge) 

At this location, the flow was investigated at the leading edge. The measured 
velocity is illustrated in Fig. B.2 as a function of the plunge phase angle, $. 

3. Third Location, X/C = 0.25 

The measured velocity is illustrated in Fig. B.3 as a function of the plunge phase 
angle, cf>. 

4. Fourth Location, X/C = 0.5 

The measured velocity is illustrated in Fig. B.4 as a function of the plunge phase 
angle, cf>. 

5. Fifth Location, X/C = 0.75 

The measured velocity is illustrated in Fig. B.5 as a function of the plunge phase 
angle, (f>. 
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Figure B.l. Measured Axial Velocity, u, at X/C=-0.25 
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Figure B.2. Measured Axial Velocity, u, at the leading edge 
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Figure B.3. Measured Axial Velocity, u, at X/C=0.25 
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Figure B.4. Measured Axial Velocity, u, at X/C = 0.5 
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Figure B.5. Measured Axial Velocity, u, at X/C =0.75 
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Plunge Phase Angle, 0, [deg] 


Z/C= -1.0063 
Z/C=-0.7813 
Z/C= 0.81165 



Plunge Phase Angle, 0, [deg] 


Z/C= -1.5063 
Z/C=-1.2563 
Z/C= 1.362475 



Plunge Phase Angle, 0, [deg] 


Z/C= -0.5563 
Z/C= -0.31167: 
Z/C=0 



Plunge Phase Angle, 0, [deg] 


Figure B.6. Measured Axial Velocity, u, at X/C = 1.0 


6. Sixth Location, X/C = 1.0, (trailing edge) 

The measured velocity is illustrated in Fig. B.6 as a function of the plunge phase 
angle, qb. 

7. Seventh Location, X/C = 1.5 

The measured velocity is illustrated in Fig. B.7 as a function of the plunge phase 
angle, qb. 
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Figure B.7. Measured Axial Velocity, u, at X/C 


1.5 
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Figure B.8. Measured Axial Velocity, u, at X/C=2.0 


8. Eighth Location, X/C = 2.0 

The measured velocity is illustrated in Fig. B.8 as a function of the plunge phase 
angle, qb. 

9. Ninth Location, X/C = 3.0 

The measured velocity is illustrated in Fig. B.9 as a function of the plunge phase 
angle, qb. 
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Z/C= -2.0063 
Z/C=-1.7563 
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Z/C= -0.5563 
Z/C= -0.31167^ 
Z/C=0 



Figure B.9. Measured Axial Velocity, u, at X/C=3.0 
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APPENDIX C: UNSTEADY MEASUREMENTS 
FOR SINGLE FLAPPING WING 


155 




X/C = -0.2 


Z/C=1 
Z/C=1.2 



Plunge Phase angle 


Figure C.l. Single Flapping Wing at X/C=-0.2, k=l, u=2[m/s], f=5[Hz], Re=8760 
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x/c = o 



Z/C=-1.4 
Z/C=-1.2 
Z/C=-1 


x/c = o 



Z/C=-0.2 

Z/C=0 

Z/C=0.2 


x/c = o 



Z/C=-0.8 

Z/C=-0.6 

Z/C=-0.4 


x/c = o 



Z/C=0.4 

Z/C=0.6 

Z/C=0.8 


x/c = o 


Z/C=1 
Z/C=1.2 
Z/C=1.4 



Figure C.2. Single Flapping Wing at X/C=0, k=l, u=2[m/s], f=5[Hz], Re=8760 
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X/C = 0.2 
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Plunge Phase angle 


X/C = 0.2 



Z/C=-0.2 

Z/C=0 

Z/C=0.2 




X/C = 0.2 


Z/C=0.4 

Z/C=0.6 

Z/C=0.8 




Plunge Phase angle 


Figure C.3. Single Flapping Wing at X/C=0.2, k=l, u=2[m/s], f=5[Hz], Re=8760 
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X/C = 0.6 



X/C = 0.6 



X/C = 0.6 



Plunge Phase angle 


X/C = 0.6 



X/C = 0.6 



Z/C=1 
Z/C=1.2 
Z/C=1.4 



Figure C.4. Single Flapping Wing at X/C=0.6, k=l, u=2[m/s], f=5[Hz], Re=8760 
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X/C = 1 



X/C = 1 



X/C = 1 



X/C = 1 



X/C = 1 


Z/C=1 
Z/C=1.2 
Z/C=1.4 



Figure C.5. Single Flapping Wing at X/C=1.0, k=l, u=2[m/s], f=5[Hz], Re=8760 
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APPENDIX D: UNSTEADY MEASUREMENTS 
FOR THE FLAPPING WING IN PRESENCE 
OF GROUND PLANE 
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X/C = -0.25 


u/U at Z/C= -0.79 
u/U at Z/C=-0.69 



Plunge phase angle 


X/C = -0.25 


u/U at Z/C=-0.25 
u/U at Z/C=0 



Plunge phase angle 


X/C = -0.25 


u/U at Z/C=0.8 
u/U atZ/C=1.3 



Figure D.l. Flapping Wing in Ground Effect at X/C=-0.25 


162 


















x/c = o 


u/U at Z/C= -0.79 
u/U at Z/C=-0.69 



Plunge phase angle 


u/U at Z/C=-0.25 
u/U at 27C=0 



Plunge phase angle 


X/c — o ♦ u/U at 27C=0.3 

-■- u/U at Z/C=0.8 

-^- u/U at Z/C=1.3 



Figure D.2. Flapping Wing in Ground Effect at X/C=0 
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X/C=0.25 - — m - u/U at Z/C=-0.79 

- <*■ , . u/U atZ/C=-0.69 



Plunge phase angle 


X/C=0.25 -□- u/U atZ/C=-0.25 

u /UatZ/C=0 



Plunge phase angle 


X/C=0.25 -■- u/U at Z/C=0.3 

-*- u/U at Z/C=0.8 



Plunge phase angle 


Figure D.3. Flapping Wing in Ground Effect at X/C=0.25 


164 







































Plunge phase angle 



Plunge phase angle 


X/C = 0.5 



u/U at Z/C=0.3 
u/U at Z/C=0.8 
u/U at Z/C=1.3 



Figure D.4. Flapping Wing in Ground Effect at X/C=0.5 
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X/C = 0.75 


u/U at Z/C= -0.79 
u/U at Z/C=-0.69 
u/U at Z/C=-0.59 



Plunge phase angle 


X/C = 0.75 


u/U at Z/C=-0.25 
u/U at Z/C=0 
u/U at Z/C=0.25 



X/C = 0.75 


u/U at Z/C=0.3 
u/U at Z/C=0.8 
u/U at Z/C=1.3 



Figure D.5. Flapping Wing in Ground Effect at X/C=0.75 
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Plunge phase angle 
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u/U at Z/C=-0.25 
u/U at Z/C=0 



Plunge phase angle 


X/C = 1 


u/U at Z/C=0.3 
u/U at Z/C=0.8 



Plunge phase angle 


Figure D.6. Flapping Wing in Ground Effect at X/C=1.0 


167 































X/C = 1.5 


u/U at Z/C= -0.79 
u/U at Z/C=-0.69 
u/U at Z/C=-0.59 



Plunge Phase Angle 


X/C = 1.5 


u/U at Z/C=-0.25 
u/U at Z/C=0 
u/U at Z/C=0.25 



X/C = 1.5 


u/U at Z/C=0.3 
u/U at Z/C=0.8 
u/U at Z/C=1.3 



Figure D.7. Flapping Wing in Ground Effect at X/C=1.5 
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- « - u/U at Z/C= -0.79 

X/C = 2 —»-u/U at Z/C=-0.69 

-• -i r - u/U at Z/C=-0.59 



Plunge phase angle 



—•— u/U at z/C-0.3 

X/C = 2 -□- u/U at Z/C=0.8 

-*- u/U at Z/C=1.3 



Figure D.8. Flapping Wing in Ground Effect at X/C=2 
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